Superinsulator as a phase of bi-particle localized states 
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We propose a physical picture of superinsulator observed recently in experiments with super- 
conducting films in a magnetic field. On the basis of previous numerical studies we argue that 
a moderate attraction creates bi-particle localized states at intermediate disorder strength when 
noninteracting electron states are delocalized and metallic. Our present numerical study show that 
such localized pairs are broken by a static electric field which strength is above a certain threshold. 
We argue that such a breaking of localized pairs by a static field is at the origin of superinsulator 
breaking with a current jump observed experimentally above a certain critical voltage. 

PACS numbers: 74.20.-z, 71.30.+h, 72.15.Rn 



I. INTRODUCTION 



An interplay of disorder, Anderson localization and su- 
perconductivity attracts active experimental and theo- 
retical interest (see e.g. reviews [1, 2] and a recent re- 
search article [3]). A weak disorder does not significantly 
affect the superconducting phase in agreement with the 
Anderson theorem [4, 5]. However, a relatively strong 
disorder can lead to a nontrivial situation when an at- 
traction between electrons creates bi-particle localized 
states (BLS phase) from noninteracting metallic delocal- 
ized electron states [6, 7] (see Fig. la,b). First signa- 
tures of this BLS phase has been obtained in the frame 
of the generalized Cooper problem [8] of two interact- 
ing particles above a frozen Fermi sea in presence of dis- 
order. Further extensive quantum Monte Carlo studies 
of two and three-dimensional (2D, 3D) Hubbard model 
[9, 10], with attraction, disorder and about hundred elec- 
trons, confirmed the attraction induced picture of local- 
ized Cooper pairs proposed in [6, 7]. This picture is in 
a qualitative agreement with 2D disordered films exper- 
iments on superconductor-insulator transition (SIT) of 
Gantmakhcr et al. [11], which clearly display a presence 
of localized pairs in the insulating phase appearing from 
the superconducting phase at relatively large magnetic 
field (see also [2] and Refs. therein). At even larger mag- 
netic fields these samples show a metallic type behavior 
which is argued to correspond to underlying metallic non- 
interacting states [7]. 

Recent experiments also discovered that in 2D disor- 
dered films the above insulating phase is abruptly de- 
stroyed by a static dc-voltage [12-16]. It was argued 
that this unusual insulating phase is related to a certain 
collective state named superinsulator [15]. The physi- 
cal origin of this state is under hot theoretical debates 
[15, 17-20]. However, this physical problem involves non- 
trivial interplay of interactions, disorder and localization 
which is not easy to handle by purely analytical methods. 
Indeed, it is known that repulsive or attractive interac- 
tions can produce derealization of two interacting par- 
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FIG. 1: Derealization of two interacting particles with at- 
tractive Hubbard interaction U = —2V by a static electric 
force F. The generalized Cooper problem is considered on 
2D lattice of size L x L = 40 x 40 at disorder strength 
W — 5V, a static field F is directed along y— axis. Proba- 
bility is shown for a lowest energy eigenstate with a maxi- 
mum probability f(y) = £^ f(x,y) at y — L/2. Left panels 
show one-particle probability f(x, y) and right panels show 
interparticle distance probability fd(x,y). The static electric 
force, directed along y— axis, is F = 0(a, b), F = 0.003L(c, d), 
F = 0.016V(e, /), F = 0.052L(g, h). 



tides above Fermi level when all one-particle states are 
exponentially localized due to the Anderson localization 
(see e.g. [21-26]). This two interacting particles (TIP) 
effect leads to an effective 3D Anderson transition for 
TIP excitations at certain energy above the Fermi level 
in the case of Coulomb or other long range interactions 
[27, 28]. But at the same time attractive interactions 



2 



create the BLS phase in the vicinity of Fermi level even 
when noninteracting states are metallic and delocalized 
[6, 7, 9, 10]. Due to that for a better understanding 
of physics of superinsulator it would be rather useful to 
study numerically the effects of a static electric held on 
localization-delocalization transition in presence of inter- 
actions. A certain progress has been reached in studies 
of attractive Hubbard model with disorder by quantum 
Monte Carlo methods since there is no sign problem in 
such a case and numerical simulations can be done with 
a large number of electrons (see e.g. [9, 10, 29] and Refs. 
therein). However, this approach is not easy to adapt to 
a case of finite static field. 

In this work we study effects of static field in the frame 
of the generalized Cooper problem using the approach 
of two interacting particles with an attractive Hubbard 
interaction U [6, 7]. Qualitative description of a physi- 
cal picture of BLS phase and its destruction by a static 
electric field F is given in Section II. In Section III we 
describe our model of the generalized Cooper problem in 
presence of rfc-ficld, the results of numerical simulations 
are presented in Section IV, discussion of physical results 
and comparison with experiments are given in Section V. 

II. PHYSICAL PICTURE OF SUPERINSULATOR 

The results obtained in [6, 7, 9, 10] give the following 
qualitative physical picture of the BLS phase, which is at 
the origin of superinsulator as it is argued below: 

i) For a moderate disorder, which e.g. by a factor two 
or less smaller than the critical value for the Anderson 
transition for 3D case, noninteracting electron states are 
still delocalized and metallic. A similar situation appears 
for finite size 2D samples where one-particle localization 
length is larger than the sample size. 

ii) However, an attractive Hubbard interaction creates 
singlet spin pairs of two electrons which total mass is 
twice larger than electron mass, hence, an effective hop- 
ping matrix element of such a pair becomes twice smaller 
or, equivalently, effective strength of disorder becomes 
twice larger that leads to localization of pairs and BLS 
phase. Such localized pairs have a certain localization 
size £ and a coupling energy A. It is important to stress 
that the Bogoliubov - de Gennes meanfield approach is 
not able to capture this BLS phase [6, 7, 9, 10]. The BLS 
phase is an insulator existing at temperature T < A. In 
a certain sense attraction between two electrons creates 
an effective well which captures localized electron pairs. 
However, excitations with energy AE > Ed ~ A are 
delocalized. Indeed, even if all one-particle states are 
localized, the excited states above the Fermi level be- 
come delocalized by interactions between electrons as it is 
shown in [24, 26-28]. Moreover, the energy excitations of 
BLS phase are delocalized above certain energy Ed ~ A 
since noninteracting states are metallic. According to the 
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FIG. 2: (Color online) Phase diagram in the 3D Anderson 
model with disorder strength W and Hubbard attraction U: 
superconductor (S), phase of localized pairs (BLS), Anderson 
insulator (7); W c ~ 16. 5V is the critical disorder for nonin- 
teracting particles at half filling. The data, taken from [6], 
are obtained in the frame of the generalized Cooper problem 
with TIP. The dashed curve with arrow schematically shows 
variation of effective U and W values with the increase of 
magnetic field corresponding to the experiment [11] (see Sec- 
tion II for more details). This phase diagram is shown for 
zero temperature T. 

Fermi-Dirac statistical distribution an excitation proba- 
bility to high energy drops exponentially and therefore 
for T < A the resistivity R xx is characterized by the 
Arrhenius law R xx oc exp(To/T) with To ~ Ed ~ A. 

Hi) In real 2D superconducting films, studied experi- 
mentally (see e.g. [2, 11]), an increase of magnetic field 
B up to a few Tesla effectively decreases attraction in- 
side Cooper pairs due to breaking of time reversal and 
also effectively increases an effective strength of the dis- 
order since magnetic field increases a return probability 
and scattering on impurities. As a result superconduc- 
tivity disappears, electron pairs become localized, but 
at even stronger magnetic field attraction between elec- 
trons is completely eliminated and one obtains a metal 
of noninteracting electrons. This effective change of at- 
traction U and disorder strength W with an increase of 
magnetic field B is schematically shown in Fig. 2 by a 
dashed curve. In this picture, resistivity R xx initially 
grows with increase of B but above a certain value it 
starts to decrease with B giving a peak in R xx , which is 
a characteristic feature of experiments (see e.g. Fig. 1 in 
[11])- 

iv) Thus, the BLS phase and its energy excited states 
can be viewed as a sequence of wells with localized pairs 
and continuum of delocalized states as it is schematically 
shown in Fig. 3a. A static electric field with force F 
creates an energy slope leading to a tilted washboard po- 
tential as it is schematically shown in Fig. 3b. Above 
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FIG. 3: (Color online) Pictorial image of the energy spectrum 
of BLS pairs localized by a Hubbard attraction and located 
at various positions along y~ direction of the lattice at static 
force F = (a); derealization of BLS pairs by a washboard 
potential V w = —Fy of static field at F > F c (b). 

a certain critical force F > F c the localized states in- 
side a well (localized pairs) become coupled to continuum 
states of dclocalizcd electrons that creates an avalanche 
of delocalized electrons. For F > F c all electrons be- 
come delocalized that produces a sharp increase (jump) 
of electron current which is a characteristic feature of su- 
perinsulator experiments [12-16]. The critical field F c 
above which the superinsulator is destroyed can be esti- 
mated by taking into account that the coupling energy of 
localized pairs A should be comparable with the energy 
change in a static field F c on the pair size i that gives 

F c nA/l, V C = F C L (1) 

The physical meaning of this relation is rather direct: a 
strong field breaks BLS pairs and creates a charge cur- 
rent. A critical voltage V c , at which a jump of current 
takes place in an experiment, is proportional to the sam- 
ple size L. 

In next Sections we present numerical simulations of 
the generalized Cooper problem in a tilted potential 
which justifies this physical picture. 



III. MODEL DESCRIPTION 

For our numerical studied of delocalization of BLS 
pairs by a static electric field F we use 2D Anderson 
model with disorder strength W and Hubbard attraction 
U between two particles. Following [7], with the same 



notations, we use the one-particle Hamiltonian: 

#i = £(^ + F - n )l n >H + ^ E l n >< n 'l ( 2 ) 

n (n,n'> 

where n and n' are index vectors on the two- 
dimensional square lattice with periodic boundary con- 
ditions in ir— direction, and zero boundary conditions in 
y— direction, V is the nearest neighbor hopping term and 
the random on-site energies E n are homogeneously dis- 
tributed in the energy interval [— ^1 . We choose the 
direction of a static electric force F along y— axis. We 
consider a square lattice with linear size L up to L = 40. 
At such sizes the eigenstates at a half filling v = 1/2 are 
practically delocalized over the whole lattice for W < 7V 
and F = so that such finite samples can be considered 
to be metallic (see more details in [7]). We consider the 
particles in the singlet state with zero total spin so that 
the spatial wavefunction is symmetric with respect to 
particle permutation (interaction is absent in the triplet 
state). 

To take into account the effects of Hubbard interaction 
we write the TIP Hamiltonian in the basis of noninter- 
acting eigenstates at F = 0: 

m i +E m2 )x ffli ,7712 

m 1 

F„ l2t m' Xmi,m') + U E Qmi^iV^m^m'^mj 
m 1 ,m 2 

= EXmi,m 2 - (3) 

Here E m are the one-particle eigenenergies correspond- 
ing to the one-particle eigenstates \4>m) and Xm 1 ,m 2 
are the components of the TIP eigenstate in the non- 
interacting eigenbasis |</> mi , 4>m 2 ) at F — 0. The matrix 
elements F mi ^ m i describe the static force transitions be- 
tween one-particle eigenstates \4> mi , <t>m 2 )' The matrix el- 
ements UQ mi m2 m > give the interaction induced tran- 
sitions between non- interactive eigenstates \<j> mi > 4>m 2 ) 
and \(f> i ,<f) i ). These matrix elements are obtained by 
rewriting the Hubbard interaction in the non-interactive 
eigenbasis of model (2) at F = 0. In the analogy 
with the original Cooper problem [8] the summation in 
(3) is done over the states above the Fermi level with 
eigenenergies E m > > Ep with m 1 2 > 0. The Fermi 

energy Ep m is determined by a fixed filling factor 
v = 1/2. To keep the similarity with the Cooper prob- 
lem we restrict the summation on m 1 2 by the condition 
1 < m 1 + m 2 < M . In this way the cut-off with M un- 
perturbed orbitals introduces an effective phonon energy 
E ph = huj D « 3.75VM/L 2 = 3.7W/a where L is the 
linear system size. When varying L we keep a — L 2 /M 
fixed so that the phonon energy is independent of sys- 
tem size. All the data in this work are obtained with 
a = 15 but we also checked that the results are not sen- 
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FIG. 4: Dependence of the average IPR £ of lowest energy 
states on the static force F; average is done over Ni = 15 
lowest energy states with a maximum of f(x, y) inside a stripe 
0.4L < y < 0.6L in a middle of the lattice at y — L/2 for 
Nd = 30 disorder realizations; the lattice size is L = 20 
(•), L — 30 (A), L = 40 (■) at U = -IV; the values of £ 
for noninteracting case are shown by open symbols (□) for 
L = 40, U = 0. Here W = 5V. 



sitive to the change of a. We note that the Hamiltonian 
(3) exactly describes the noninteracting problem. 

To analyze the effects of static force on localized 
BLS pairs we solved numerically the Schrodingcr equa- 
tion (3). After that we rewrite the obtained eigen- 
states in the original lattice basis with the help of the 
relation between lattice basis and one-particle eigen- 
states |n) = J2 m Rn, m \<p m ) ■ As a result of this pro- 
cedure we obtain the two-particle probability distribu- 
tion /2 (111,112) from which we extract the one-particle 
probability /(n) = J2n 2 h( n i, n 2) and the probability 
of interparticle distance /<j(r) = En 2 / 2 ( r + n2 ' n2 ) w ^ n 
r = ni — n.2- The localization properties are character- 
ized by the one-particle inverse participation ratio (IPR) 

£ = E„/(n)/E„/ 2 (n)- 

While in [7] only the ground state properties of given 
disorder realization have been studied here we also in- 
vestigate the properties of excited states with the TIP 
energy AE counted from the Fermi level of noninteract- 
ing particles: AE = E — 2Ef, where E is the eigenen- 
ergy of (3). We also consider only eigenstates with a 
maximum of one-particle probability inside the space 
range — L/4 < y < L/4 to avoid finite size effects in 
y— direction. In addition, we analyze only those eigen- 
states where the overlap probability of TIP to be on the 
same site is relatively large /d(0, 0) > 5/L(2L— 1). In this 
way the states with strongly separated particles are elimi- 
nated. Such an approach approximately corresponds to a 



finite particle density. We use usually Njj = 30 disorder 
realizations for statistical average. 

Below we present numerical results for U = — 2V, 
W = 5V and L < 40 at various values of F. The de- 
tailed studies presented in [7] ensure that at F = these 
conditions are located well inside the BLS phase when 
the noninteracting states are delocalized (see Fig. lb in 
[7]) while the ground state in presence of Hubbard at- 
traction is well localized (see Fig. le,f in [7] and Fig. la,b 
here). We checked that the behavior in F remains similar 
at other values of parameters, e.g. W = 3V. 

IV. NUMERICAL RESULTS 

The derealization of BLS pairs by a static electric 
force F is illustrated in Fig. 1 for one specific disorder 
realization. Here the disorder strength is relatively weak 
so that at U — noninteracting eigenstates taken at half 
filling v = 1/2 and Ep rts are delocalized over the whole 
lattice of size L = 40 (see Fig. lb in [7]). At F = a 
moderate Hubbard attraction U = —2V creates local- 
ized pairs with a certain coupling energy A (Fig. la,b). 
This localization remains robust against a weak static 
field (Fig. lc,d) but at larger fields the localization is de- 
stroyed by a static force and particles become delocalized 
over the whole lattice (Fig. le,f,g,h; to avoid boundary 
effects at finite F we select states in the middle of the 
lattice at y = L/2). The probability to have two parti- 
cles close to each other also drops drastically for F > F c 
clearly demonstration pair breaking. 

The dependence of average IPR £ of lowest energy 
states on the static force F is shown in Fig. 4. At small 
F < F c the values of £ are size independent being much 
smaller compared to the case of noninteracting particles. 
This shows that a Hubbard attraction gives localization 
of pairs at low energy. For F > F c IPR starts to grow 
with the system size L demonstrating breaking of pairs 
and particle derealization over the whole lattice. Ac- 
cording to the data of Fig. 4 we have F c « 0.015V at 
given U = —2V and W = 5V. The data in Fig. 5 give 
the coupling energy A w 0.1V. Hence, this F c value 
is in a good agreement with a simple estimate (1) with 
a numerical factor A = F c £/A w 1 corresponding to 
£{F = 0) » y/£(F = 0) w 6 and A = 0.1V. 

The dependence of £ on coupling energy AE = E — 
2Ep, for states in the middle of the lattice at y w L/2, 
is shown in Fig. 5. According to this data the coupling 
energy is A w 0.1V. This value is by a factor 2 smaller 
than the one found in [7] since only one lowest state for a 
given disorder realization was taken in [7] while here we 
average over few lowest states and also allow a relatively 
weak overlap fd between TIP states. The states with 
— A < AE < are well localized at F — since its 
IPR £ is independent of lattice size L. In contrast, for 
F = 0.036F > F c the IPR £ grows with the system size 
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FIG. 5: (Color online) Average IPR £, for states peaked inside 
a stripe OAL < y < 0.6L , as a function of the coupling energy 
AF = E - 2E F for U = -2V, W = 5V and lattice size L = 20 
(blue circles), L = 30 (red triangles), L — 40 (black squares), 
at field F = (open symbols) and F = 0.036V (full symbols). 
The same averaged IPR £ for 17 = 0, W = 5V, F = and 
L = 40 is shown by (+) symbols. 




FIG. 6: (Color online) Average IPR £, for states peaked inside 
a stripe 0.4L < y < 0.6L , as a function of the coupling energy 
AE = F-2F F and static force F for (7 = -2V, IV = 5V (top 
panel) and (7 = 0, W = 5V (bottom panel); here L = 40. 



L showing that in this regime the states are delocalized. 

At energies AE > the IPR grows with energy and be- 
comes comparable with the system size L 2 ; also its is not 
sensitive to F. This is in agreement with the fact that 
noninteracting states are delocalized. Also interaction 
for excited states gives an additional TIP derealization. 
For large values of L — 40 the static force gives a cer- 
tain restriction of eigenstates spreading along the force 
direction due to the TIP energy conservation that gives 
a decrease of IPR value comparing to the case F = 0. 

A more detailed dependence of £ on coupling energy 
AE and static force F is given in Fig. 6. For U — we 
have £ « 200 which is practically independent of AE and 
F while for U = — 2V we have very small £ ~ 10 for F = 
and large £ - 300 for large F at -0.1 < AE/V < 0. 
For energies AE > the states are delocalized at all 
F. These data also confirm the picture of field induced 
destruction of the BLS phase and derealization. 

Dependence of £ on y and coupling energy AE is shown 
in Fig. 7 for different values of F with and without inter- 
action. For U = we have £ practically independent of y 
and AE in agreement with previous data of Figs. 4,5,6. 
In contrast, in presence of attraction the states in the 
middle of the lattice in y have small £ (localized) at low- 
est energies for small F (panels b,d) and have large £ (de- 
localized) for large F (panels f,h). However, at the ends 
of the lattice in y direction the values of £ are less sensi- 
tive to F due to boundary effects. Indeed, the static field 
forces particles to stay close to the boundary at y— ends 
of the lattice and hence the field induced derealization 
is not well visible in this region. Due to that reason we 
use the states in the middle of the lattice to detect field 
induced derealization in a clear way in Figs. 4,5,6. 



V. DISCUSSION 

The obtained numerical results confirm the physical 
picture of superinsulator destruction by a static field de- 
scribed in Section II: the BLS pairs, localized by attrac- 
tion inside noninteracting metallic phase, in presence of 
static field start to be coupled with higher energy delo- 
calized states and above certain threshold F > F c (1) all 
pairs become delocalized. This creates an avalanche of 
delocalized electrons which gives an enormous increase of 
current in agreement with experimental observations. Of 
course, our numerical data detect derealization of only 
one pair in a given disorder realization. However, the dis- 
tribution of values of pair coupling energy A = 2Ep — E g 
is strongly peaked near its average value (see Fig. 8) so 
that a large fraction of localized pairs becomes delocal- 
ized approximately at the same static field that gives a 
sharp current growth for V > V c = F C L (here E g is a 
ground state energy for a given disorder realization at 
F = 0). 

At that point we would like to note that even if our at- 
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FIG. 7: (Color online) The IPR £ for W = 5V at U = 
(left column) and t/ = — 2 V (right column), at different 
values of the static force F = (a,b), F = 0.004V (c,d), 
F = 0.016V (e,f), F = 0.04V (g,h). Each panel has 40 x 40 
cells, the vertical direction corresponds to the coupling en- 
ergy AE = E — 2Ef, the horizontal one to 1 < y < 40. 
The bottom row corresponds to the lowest coupling energy, 
and the upper row corresponds to the highest coupling en- 
ergy within the energy intervals at U — being (0, 0.54V) (a), 
(-0.05V, 0.59V) (c), (-0.25V, 0.78V) (e), (-0.69 V, 1.18V) (g) 
and at U = -2V being (-0.2V, 0.39V) (b), (-0.23V, 0.57V) 
(d), (-0.39V, 0.75V) (f), (-0.77V, 1.17V) (h). The cell color 
gives the average £ inside the cell (with AE being in the corre- 
sponding energy range defined by the row, and the maximum 
of probability distribution f(y) along y, being located at y po- 
sition defined by the column). For each panel ~ 30000 states 
are used with Na = 30 disorder realizations. These states are 
selected in such a way that the probability of two particles 
located at the same site is greater than 5/L(2L — 1); here 
L = 40. 
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FIG. 8: (Color online) Distribution histogram of pair coupling 
energy A = 2Ef — E g obtained from Nd = 1000 disorder 
realizations at U = -2V, W = 5V, F = 0, L = 40; here N 
gives a number of realizations found in a given cell of A/V, 
Eg is the ground state energy of a given realization obtained 
numerically from (3). 



traction is formally relatively strong (e.g. \U\ = 2V ~ V) 
it effectively gives a rather weak coupling energy of BLS 
pairs A « 0.1V (see Fig. 8). There are a few physical 
reasons behind this. At first, a simple physical estimate 
gives a consistent value A ~ \U\/£ ~ 0.05V <C V where 
we use numerically found value of IPR £ « £ 2 ~ 40 from 
Fig. 4. Thus the numerical value of A is by a factor 100 
smaller than the energy band width of the noninteract- 
ing 2D problem E b w 8V + W w 14V. Second, in (3) 
the summation over one-particle orbitals is done only in- 
side a Debye energy interval E p h = Hujd ~ 3.75V/a w 
0.25V -C Eb and due to that the effective attraction is 
additionally decreased giving a relatively small coupling 
energy. Also in a good metallic phase (e.g. W < 2V) 
the interaction (e.g. even \U\ = 2V) produces a quite 
weak effect on noninteracting delocalized states accord- 
ing to numerical results presented in [6, 7] and according 
to usual theoretical estimates for interaction matrix ele- 
ments U s ~ U / g <C Eb, which are inversely proportional 
to a sample conductance g ^> 1 (see e.g. [22, 24] and 
Refs. therein). Due to these reasons we think that the 
claim expressed in [3] , that the BLS pairs appear only as a 
result of nonphysically strong attraction, is not justified 
since our BLS pairs have rather small coupling energy 
A -c Eb and have rather large size £ 1. In addi- 
tion, the numerical data presented in [6] (see also Fig. 2) 
show that the BLS pairs appear also at smaller values of 
attraction. Thus, even if, with the aim to have well lo- 
calized pairs inside a system size available for numerical 
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FIG. 9: (Color online) Critical field F c for superinsulator de- 
struction as a function of sample size L from experiments [14] 
(points from Fig.5 there). The theory (1) is shown by red 
horizontal line. 



simulations, we fixed an attraction at a relatively strong 
value, our numerical data show that such a choice is still 
in the regime of weakly coupled pairs of relatively larger 
size that corresponds to the experimental regime. We 
note that even larger |£/|/V values are typically used in 
quantum Monte Carlo simulations (see e.g. [29]). 

The main result of this studies is given by Eq. (1) which 
determines the critical voltage V c of superinsulator de- 
struction via the sample size L, BLS coupling energy A 
and pair size I. According to (1) the critical voltage 
V c is proportional to the sample size and hence, F c is 
independent of L. This is in a good agreement with the 
experimental data obtained in [14] (see Fig. 9 with exper- 
imental points from Fig.5 in [14]). Indeed, in the range 
0.5/im < L < 150/im the value of F c shows significant 
fluctuations but in average remains constant. The aver- 
age value is F c ss 30V/ cm, and since the typical value of 
pair coupling energy is A ks Tq ~ 3 — 15K, we find the 
size of localized pairs to be £ = Tq/F c = 100 — 500tito. 
The theory (1) is valid for L > I where indeed F c is 
independent of L, a part of fluctuations. However, for 
L < £ ~ 0.5/im one enters into another regime where 
the sample size becomes comparable with the pair size 
that can lead to an increase of F c seen experimentally. 
It is clear that as soon as the superinsulator phase is de- 
stroyed at F > F c a further decrease of F below F c places 
localized pairs in other new locations where due to fluc- 
tuations of disorder one gets a comparable but somewhat 
different new value of F c . This leads to a hysteresis be- 
havior observed experimentally. 

It is also interesting to note that the relation (1) allows 
to determine the dependence of i on magnetic field. In- 
deed, in experiments [12, 13, 15] the values of L is known, 
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FIG. 10: (Color online) Schematic phase diagram in the 
temperature-disorder plane (T, W) for quasi-2D or 3D Ander- 
son model at fixed Hubbard attraction U = —2V and fixed 
filling factor with the Fermi energy Ef- the vertical dashed 
line shows the Anderson transition at W = W c ~ 16.5V 
for noninteracting particles, S is the superconducting phase, 
BLS is the insulator phase of localized BLS pairs, M is essen- 
tially the noninteracting metallic phase, TIP is the metallic 
phase of delocalized TIP pairs. Here we assume T <C Ef, 
static field is zero, temperature linear axis is shown in arbi- 
trary units. 



and also A w To and V c are experimentally known as a 
function of magnetic field B that allows to determine the 
dependence £(B) from the relation (1). 

On the basis of presented results we can draw a global 
phase diagram in the temperature-disorder plane shown 
in Fig. 10 for a fixed attraction strength (e.g. at U = 
—2V), zero static field and fixed Fermi energy Ep. At 
small disorder and temperature we have the supercon- 
ducting phase S which is followed by a transition to 
metallic phase M at large temperature or to the local- 
ized BLS phase at low temperature. At a larger disorder 
W > W c , but still small temperatures, the BLS phase en- 
ters in the insulating regime of noninteracting Anderson 
insulator I. In this regime with W > W c but tempera- 
ture above a certain threshold T > T 2 ~ |f7|(l/£ + l/£i), 
the TIP pairs become delocalized by interactions with 
emergence of metallic phase of TIP delocalized pairs, as 
it is discussed in [21, 22, 24]. Here £i is a noninteracting 
one-particle IPR which gives a dominant contribution for 
a disorder W > W c which is not very close to the critical 
point W = W c , the term l/£ with IPR of BLS is included 
to have interpolation between two phases. 

In conclusion, we presented the BLS based physical 
picture for a destruction of superinsulator by a finite 
static field observed experimentally in [12-16]. This pic- 



8 



ture is rather different from other theoretical explana- 
tions discussed in [15, 17-20]. The main new element 
of our theory is the existence of localized pairs created 
by attraction inside noninteracting metallic phase which 
is absent in the above theoretical proposals. In our the- 
ory the BLS phase is the basis of superinsulator and since 
the noninterating states are metallic this phase is sharply 
broken by a finite static field which breaks electron pairs, 
localized by attraction, and lets them propagate like non- 
interacting particles in a metal. Our analysis considers 
a breaking of only two interacting particles. In a real 
system with many pairs a breaking of a few pairs can 
create a strong avalanche and breaking of other pairs so 
that a critical static field can be determined by break- 
ing of mostly weakly coupled pairs with a smaller crit- 
ical fields, compared to average field values found here. 
This physical picture is rather different from the superin- 
sulator picture discussed in [15]. However, we keep the 
term superinsulator which in our opinion nicely describes 
impressive experiments [12, 13, 15, 16] with supercon- 
ducting films, which become insulating in magnetic fields 
with abrupt emergence of current above a certain critical 
dc— voltage. 
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